function xx=fun2(t,x) 
r1=1;k1=0.1;r2=0.5;k2=0.02;
global e;
xx=[x(1)*(r1-e-k1*x(2));x(2)*(-(r2+e)+k2*x(1))];
